Gapless edge states and their stability in two-dimensional quantum magnets 
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We study the nature of edge states in extrinsically and spontaneously dimerized states of two- 
dimensional spin-i antiferromagnets, by performing quantum Monte Carlo simulation. We show 
that a gapless edge mode emerges in the wide region of the dimerized phases, and the critical 
exponent of spin correlators along the edge deviates from the value of Tomonaga-Luttinger liquid 
(TLL) universality. We also demonstrate that the gapless nature at edges is stable against several 
perturbations such as external magnetic field, easy-plane XXZ anisotropy, Dzyaloshinskii-Moriya 
interaction, and further-neighbor exchange interactions. The edge states exhibit non TLL behavior, 
depending strongly on model parameters and perturbations. Possible ways of detecting these edge 
states are discussed. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Kt, 73.43.-f, 75.10.-b 
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Introduction.— In recent years, gapful ground states 
without any local order parameter and their boundary 
properties have been vividly studied in quantum many- 
body physics from both theoretical and experimental 
viewpoints. Among such disordered states, for instance, 
two- and three-dimensional (3D) topological insulators 
(TIs) have attracted much attention as novel many-body 
phases in solids^ 2 - Their fundamental properties are that 
the bulk has a finite excitation gap, but its boundary 
(surface or edge) is metallic (i.e., gapless) and around 
the boundary, up-spin electrons move antiparallcl to the 
motion of down-spin electrons (this nature is called heli- 
cal). This gapless boundary state is quite stable against 
any perturbation with time reversal symmetry, and the 
existence of a helical edge mode is recorded in a Z 2 topo- 
logical invariant defined on the bulk (bulk-edge corre- 
spondence). 

In quantum spin systems, the Haldanc-gap stated is 
also famous as a state without any local order. This state 
is defined as the ground state of ID spin-1 antifcrromag- 
netic (AF) chains and is actually realized in some quasi- 
1D magnets4^ Its characteristic features can be captured 
by valence-bond solid (VBS) pictured Namely, the Hal- 
dane state can be approximated by the uniform tensor 
product state of local singlet dimers comprised from two 
fictitious S = j spins on neighboring sites, which are 
generated via the decomposition of original S = 1 spin 
on each site. From the solid singlet distribution, we can 
easily understand the existence of a finite excitation gap 
(called Haldane gap) on the Haldane state. Similarly, the 
uniform alignment of singlets indicates the absence of any 
local order parameter, but we can construct a non-local 
string order parameter— to distinguish the Haldane state 
from the other paramagnetic phases. The VBS picture 
also shows that an almost free S = \ spin appears at the 
edge of finite-size spin-1 Haldane state under free bound- 
ary condition.— In addition to these results, recently new 
ways of characterizing the Haldane states have been ac- 
tively discussed based on symmetries and artificial quan- 



tities such as entanglement spectra.— — 

All of the gapped, topological phases in free fermion 
systems including TI have been succeeded in being clas- 
sified theoretically*!^— On the other hand, topological 
phases and boundary states in quantum spin systems 
have been less understood, expect for a few VBS (such as 
the Haldane state) and short-range valence-bond states 
in ID spin systems including spin ladders. Therefore, 
the understanding of topological nature and boundary 
properties in quantum spin systems, especially, in higher- 
dimensional magnets, is an important, fundamental issue 
in magnetism. One might image VBS^ or exotic spin- 
liquid states as typical examples of 2D or 3D gapped non- 
magnetic spin states with a gapless edge mode. It is how- 
ever difficult to prepare such states in nature because the 
corresponding Hamiltonians contain various tuned cou- 
pling constants. Simple, realistic systems with a gapless 
edge mode would be suitable for starting to understand 
edge modes in 2D and 3D spin systems. In this paper, we 
thus study the nature of edge modes in 2D spin-Peierls 
(dimerized) states, by using quantum Mote Carlo (QMC) 
simulation based on the worm algorithm*^— 

There is a similarity between TI and the Haldane state; 
A finite bulk gap, gapless boundary states, and Z 2 topo- 
logical invariant of TI seem to correspond to a Haldane 
gap, free edge spins, and the string order parameter of 
the Haldane state, respectively. We will hence consider 
how the edge modes of dimerized states are different 
from/similar to those of TIs and the Haldane state. It 
would be impossible to define any topological order pa- 
rameter for 2D dimerized states, and in that sense, the 
gapless edge states in dimerized states (even if exist) are 
naively expected to be less stable compared to those of 
TIs. We should however note that it is generally hard 
to predict whether or not there exist gapless edge states 
and how stable they are, since quantum spin systems we 
consider below are strongly correlated systems differently 
from free fermion models for TIs. 

Model— Dimerized phases arc roughly classified into 
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FIG. 1: (color online) (a) 2D dimerized model with bond 
alternation along the y direction, (b) 2D JQs model with 
six-spin interaction term. The modified coupling constant Q' 3 
on edges is explained in the text. 



an extrinsically dimerized phase without any SSB and an 
intrinsically dimerized one with a spontaneous transla- 
tional symmetry breaking. To study those two states we 
utilize two SU(2)-symmetric spin-i AF models on square 
lattice with size LxL: The dimerized mode l 18 ' 19 and the 
JQ 3 model^ Their Hamiltonians are given as 

%dim = ^""^ JSj ■ Sj + ^""^ JaSj yy0vcn Sj ■ Sj+e y ,(l&) 

(i,3) 

(i,j,k,l,m,n) 

where Sj is the spin-^ operator on site j = (j x , j y ) (j x .y S 
Z), e x [ y ] = (1, 0) [(0, 1)] is the unit vector for the x[y] di- 
rection, J > is the AF exchange coupling constant be- 
tween neighboring spins, and CV/ = 1/4 — Si ■ Sj . In the 
dimerized model ([lap , a denotes the magnitude of exter- 
nal dimerization along the y direction as shown in Fig. [TJ 
in which the dimerized bond strength J' = J(l + a). 
If a is strong enough and the open boundary condition 
for y direction is imposed as in Fig. [TJ an effective spin 
chain is expected to appear along the x direction at the 
edge thanks to the formation of dimerization on every 
strong bonds J'. In the JQ 3 model Ijlbp. the second 
Q 3 term includes six-spin interactions, where six sites 
(i,j,k,l,m,n) are defined by two neighboring plaqucttcs 
(rectangle) shown in Fig. [TJ The symbol J2(i,j,k,l, m ,n) 
stands for the summation over all rectangles on square 
lattice. 

The ground-state phase diagrams for the dimerized 
and JQ 3 models have been investigated by QMC calcu- 
lations, and both models show the Nccl-dimcr quantum 
phase transition^ For the dimerized model, the criti- 
cal point is located at J' = J' c ~ 1.91 J 18 ' 19 and singlet 
dimers appear on all the bonds J' for J' > J' c . The tran- 
sition of the JQ3 model takes place at Q 3 = Q 3c ~ 1.5J 
and the spins spontaneously form a columnar dimer state 
along the x or y direction when Q3 > Q 3c £H It is worth 
notice that the dimer phase of Eq. (jTa)) does not ac- 
company any SSB similarly to TI, while in the case of 



Eq. (jlbl) . the translational symmetry is spontaneously 
broken. 

The JQ3 model seems to be a toy model, but it is one 
of the few spin models with a spontaneously dimerized 
phase and can be accurately analyzed by QMC simu- 
lation without negative sign problem. Furthermore, its 
dimerized ground-state wave function is expected to be 
close enough to that of real dimerized magnets. To study 
possible gapless edge modes of the model (fTa)) . we impose 
open boundary condition for the y direction as shown in 
Fig. [TJ We have confirmed from the QMC simulation 
that in the dimer phase of the JQ3 model, singlet dimers 
tend to reside on the edge. For this dimerization pattern, 
no gapless edge state is expected. In order to remove the 
dimers on edges and make the same dimerization pattern 
as that of the dimerized model (|la[) . we modify the value 
of Q 3 to Q' 3 (< Q3) at the edge as depicted in Fig. [TJ 
From QMC results of Q' 3 /Q 3 = 0, 0.25, and 0.5, we have 
checked that the resultant dimers do not reside on the 
edges and the nature of edge states is not sensitive for 
changing the value of Q' 3 . We therefore set Q' 3 /Q 3 = 0.5 
throughout this paper. 

ID liquid.— In order to judge whether or not a gapless 
edge mode are present, we utilize two-point spin corre- 
lation functions at edge of two models in Eq. ([TJ. If it 
exists, an algebraic decay of the correlators is expected, 
while the correlators decay in an exponential fashion 
when the edge state has a finite excitation gap. When 
the dimerized model approaches to the limit a — > 00, 
an isolated spin-^ AF chain appears at the edge and its 
low-energy physics is governed by a gapless Tomonaga- 
Luttinger (TL) liquid. It is therefore important to sum- 
marize spin corrrelators of the TL-liquid phase as a ref- 
erence point before embarking on our QMC results. For 
an ideal TL liquid phase of spin-i AF chains, spin cor- 
relations are known to behave as^ 

(S?SS) - Ci(-l)V-"*» + • • • , (2a) 
(S*;:^) - m 2 + C 2 r-^ cos(2fc F r) + • • • , (2b) 

at long distances r 3> 1. Here the magnetization m = 
(S') is induced by external uniform field H, 2k p = 
7r(l — 2m) is the "Fermi" wave number, and C1.2 are 
non-universal constants. It is well known that critical 
exponents r\ xy ^ z satisfy the relation r\ xy r\ z = 1^ and 
Vxy = T)z = 1 occurs at the SU(2)-symmetric case. We 
also note that a 2kp incommesurate oscillation disap- 
pears at m = H = in Eq. (|2b|) . 

Numerical analysis.— Now we discuss the QMC results 
of two models (fTaj) and (|lb|) . Figure [2] shows numerically 
determined spin correlation functions along the x direc- 
tion on edge and those along the y direction at x = L/2 
for the dimerized phases of both dimerized and JQ 3 mod- 
els at a sufficiently low temperature T. Hereafter, a cor- 
relation function along the x direction on edge (along the 
y direction) is called edge (bulk) correlation. The long- 
distance behavior of edge spin correlations S^^(r x ) = 
(S^ r 0)^(00)) are wen explained by a power-law decay 
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in both models, reflecting the existence of a gapless edge 
state. The critical exponents r\ defined by \S^(r)\ ~ r~ v 
is evaluated as n ~ 1.13 (0.97) for the dimer (JQ3) 
model. We stress that this algebraic-decay behavior sur- 
vives far from the dimer limit. In fact, the parameters 
Q3/J = 2.5 and J'/J = 2 correspond to Qs/Q3 c — 1.67 
and J'/J' c — 1.05, respectively. Furthermore, the crit- 
ical exponent 77 somewhat violates a TL-liquid relation 
Vxy = T)z — 1 • This implies the two dimensionality affects 
the gapless edge state. We have also evaluated rj ~ 1.01 
(0.99) at a deeper dimerized point J'/J = 5 (Q3/J = 4) 
in the dimerized (JQ3) model. On the other hand, as 
shown in the inset of Fig. [2] the bulk spin correlations 
B^{r y ) = (S( L/2 ry) S^ L/20) ) decay exponentially, indi- 
cating a finite dimerization gap in the bulk. The correla- 
tion lengths £, v of |B MM (r w )| ~ exp[— ry/^y] are evaluated 
as t; y ~ 7.4 (4.0) for the dimerized {JQ3) model. These 
large values of £ y clearly show that the edge chain is really 
correlated to the bulk. From these results, we conclude 
that a gapless edge state is realized in the wide range of 
both extrinsically and spontaneously spin-Pcicrls phases, 
but the exponent deviates from that of TL liquid, espe- 
cially, near the 2D critical points. 
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FIG. 3: (color online) Perturbations in the dimerized model: 
(a) DM interaction, (b) next-nearest neighbor interaction, 
and (c) bond alternation. Open boundary condition is im- 
posed and the edge state appears parallel to the x direction. 
Cross circles on J bonds indicate the direction of Di.j vector. 



The Zeeman term H. z breaks time-reversal symmetry, 
and %dm does link-parity symmetry. These two and 
Hxxz reduces the SU(2) symmetry to the axial U (1) type. 
The bond alternation term Halt eliminates the transla- 
tional symmetry along the x direction (if we apply peri- 
odic boundary condition). In contrast, Hn.n.n does not 
violate any symmetry of the original models. They are 
depicted in Fig. [3J 




FIG. 2: (color online) Spin-spin correlation ]5 , MM (r a: )] — 
I {S£ r ) S^ j ) I on the edge in the dimerized phase of (a) 
dimerized model and (b) JQ3 model. Insets are the results 
along the bulk direction: |S MM (r y )| = \{Sf L/2iry) S£ L/m )\. 

We normalize the values of |S (JM (1)| and to be unity. 

We depict only the results for B xx because of the SU (2) sym- 
metry of the models. Errors are drawn on all symbols, but 
they are less than the symbol size. 

Stability against perturbations.— We next discuss how 
robust the gapless edge mode is against various kinds 
of perturbations. This is very important since the 
lack of the stability indicates the difficulty of observ- 
ing the edge state in real magnets. As realistic per- 
turbations, we consider (1) uniform Zeeman term Ji z = 
—HJ2jSj, (2) XXZ-type magnetic anisotropy % XX z = 
- JA i{ij) S i S h ( 3 ) Dzyaloshinskii-Moriya (DM) inter- 



action %DM 



(S t x Sj) with D iS || S z , (4) 



next-nearest-neighbor interaction perpendicular to the 
edge direction Wnnn = T)i _„■ , ESi ■ Sj, and 

(5) additional bond modulation for the x direction These 
perturbations possess the following nature of symmetry. 



FIG. 4: (color online) Edge spin correlations S^^r^) of (a) 
dimerized model and (b) JQ3 model under a magnetic field 
H . Circles and triangles are the transverse (/i = x) and longi- 
tudinal (/i = z) correlators, respectively. We have eliminated 
a constant part of {Sj} 2 from the data of u = z. 
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FIG. 5: (color online) Magnetization profile m(r y ) = 
(SX/2 r ) °f the dimerized model with J'/J — 2 under mag- 
netic fields H at r x — L/2. 

First let us consider the effects of the Zeeman term. 
We plot the edge spin correlations and magnetization 
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profiles in the presence of magnetic field H in Figs. [4] 
and [5j respectively. Figure 0] shows that even in H, a 
power-low decay fashion survives in the edge spin cor- 
relations, having an incommensurate oscillation in the 
S z correlation. The stability against H is in contrast to 
the helical edge of TIs and edge spins of Haldane states. 
The incommensurability is very similar to Eq. (|2b[) , while 
such an incommensurate oscillation is absent in the bulk 
correlations. Figure [5] reveals that a finite magnetization 
emerges with increase of H only around the edge. Fi- 
nite magnetizations on sites near the edge indicates that 
some arrays around the edge also become gapless due to a 
small H. This inhomogeneous magnetization profile can 
be observed in principle, for example, by using nuclear 
magnetic resonance (NMR). Both Figs. H] and [5] clearly 
indicate that the gapless edge mode survives under field 
H. 

Next let us consider effects of 'Hdm and TL XXZ . For the 
JQs models with perturbations (2)-(5), technical difficul- 
ties of the QMC method emerge in high accurate com- 
putations. We will therefore focus only on the dimerized 
model below. However effects of the perturbations on 
the dimerized model are probably very similar to those 
on the JQ3 model since the wave functions of both extrin- 
sically and spontaneously dimerized phases are expected 
to massively overlap each other. 

Before the analysis, we should note the following prop- 
erty of the DM interaction. As we consider, for example, 
a uniform DM term with = (0,0, D z ), the Hamil- 
tonian Hdim + %dm can be mapped onto the following 
easy-plane anisotropic form 

W = ^J±(SjSj^ a + S?SV +e J + JSjSj +ex 
3 

(3) 



via a proper unitary transformation Sj: 



S 



± 23.24 



Here Jj_ = Jyl + (D z / J) 2 . Remarkably, the modified 
system ([3]) recovers the link-parity symmetry and has 
an easy-plane anisotropy. This kind of mapping can be 
applicable for a wide class of DM terms between neigh- 
boring spins if we adopt the open boundary condition for 
both x and y directions. Therefore, it is enough to study 
the dimerized model with an easy-plane XXZ anisotropy 
in order to see effects of both %dm and H xxz . 

Figure [6] (a) indicates that both the longitudinal and 
transverse edge spin correlations decay algebraically in 
the system ([3]) with D z /J = 0.1. We thus conclude 
that the gapless edge state is stable against both the DM 
and easy-plane XXZ interactions. In purely ID spin-^ 
AF chains with easy-plane anisotropy, critical exponents 
i)xy,z usually satisfy i] z > rj xy in addition to rj z j] xy = l. 22 
In the present 2D case of Fig. [S] (a), however, we obtain 
r) xy ~ 0.92 and r\ z ~ 1.48. Namely, the exponents satisfy 
the inequality rj z > r) xy , but the relation r\ z r\ xy = 1 is 
clearly broken down. The difference between critical ex- 
ponents of the ID chain and the edge of the 2D dimerized 




FIG. 6: (color online) Effects of (a) a DM interaction (equal 
to XXZ anisotropy), (b) a next-nearest neighbor interaction, 
and (c) an extra bond-alternation along the edge direction on 
the edge spin correlations of the dimerized model. Circles and 
triangles are the results of /1 = x and /j, — z, respectively. In 
panel (c), we have used semi-log scale, and have set a larger 
J' /J = 5 and a lower temperature T = 1/(2L) to depict a 
clear exponential-decay form of the edge correlation. 
JQmodel 



model would be attributed to the correlation between the 
bulk and edge. 

The edge spin correlations for the case with Hn.n.n 
are given in Fig. [6] (b), in which we have adopted a fer- 
romagnetic coupling E/J = —0.05 to avoid the negative 
sign problem. The figure shows a power-law decay of the 
correlation functions, and it indicates that the gapless 
edge mode still survives in a small Hn.n.n- We however 
note that evaluated critical exponents r] xy = r/ z ~ 0.71 
deviate from the value of ideal TL liquid similarly to the 
cases with "Hdm and H xxz . 

Figure [S] (c) is the result of the system "Hdim + "Halt- 
As expected, the edge spin correlation changes from an 
algebraic form to an exponential one due to the bond 
dimerization along the x direction. From these analy- 
ses of perturbations, we see that the gapless nature of 
the edge state survives after introducing several pertur- 
bations with different symmetries except for the bond 
alternation H a \t ■ It suggests a high possibility of the real- 
ization of a gapless edge mode in 2D spin-Peierls phases. 

Experimental ways.— Finally we discuss possible ex- 
perimental methods of probing signatures of the gapless 
edge modes in 2D spin-Peierls states. As we already men- 
tioned, a finite magnetization rapidly grows only around 
the edge sites as shown in Fig. [5j if we apply H in 
the dimerized states. Observing such a site-dependent 
magnetization (e.g., by using NMR) could be a signa- 
ture of the existence of a gapless edge state. The NMR 
relaxation rate \/T\ for nuclear spins near the edge is 
expected to contain a power-law T dependence T~ e at 
low temperatures j22 implying the existence of a gapless 
mode. Inelastic neutron scattering spectra could also 
provide information about gapless edge modes. In the 
case of TIs, ARPES spectra have been often used to see 
the gapless dispersion on the surface/edge in addition to 
gapped bulk excitations. Similarly, as shown in Fig. 
the neutron scattering spectra are expected to possess a 
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FIG. 7: (color online) Schematic inelastic-neutron-scattering 
spectrum of a 2D spin-Peierls phase with a gapless edge mode 
in the space of frequency uj and the wave number k x . The 
wave number k y is fixed to n. The weight in the low-energy 
region oj/ J < f is the contribution from the edge state, which 
lower bound would be similar to the des Cloizeaux-Peason 
mode ~ J\ sm(k x )\. The weight in the higher-energy region 
comes from the gapped triplet excitations in the bulk. 



des Cloizcaux-Peason-like gapless branc h 22 ! 25 due to the 
edge state in 2D Peierls magnets. This would be a strong 



evidence for the gapless edge mode. In addition to these 
ways, for instance, heat transport properties would cap- 
ture the nature of the gapless edge mode. 

Conclusions.— We have investigated several proper- 
ties of the edge spins of the dimerized and JQ3 models. 
When the system is in the dimerized phase and the sin- 
glet dimcrs do not reside on the edge, the gapless edge 
state really emerges. We have shown that the edge mode 
is quite robust against various perturbations with differ- 
ent symmetries. Particularly, the stability against exter- 
nal magnetic field H is in contrast with the helical edge 
modes of TIs and free edge spins of the Haldane state, 
and H induces an edge magnetization as shown in Fig. [5] 
The edge spin correlations algebraically decay like a TL 
liquid, but their critical exponents rj xy>z generally violate 
the TL-liquid relation r) xy r] z = 1, depending on the de- 
tail of the models. NMR, inelastic neutron scattering, 
and heat transport would be hopeful experimental ways 
of detecting the gapless edge modes. 

We thank J. Lou and H. Tsunetsugu for fruitful dis- 
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and cluster machines in Nano-micro structure science and 
engineering, University of Hyogo. 
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